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Abstract — The arid and semi-arid regions constitute roughly one third of the toted earth 's surface. In these regions water 
scarcity is one of the main limiting factors for economic growth. The impact of such water scarcity is amplified by inefficient 
irrigation practices, especially since about 85% of available water is used for irrigation in these regions. Therefore, a sound 
and efficient irrigation practice is an important step for achieving sustainable management of water resources in these 
regions. In this regard, a better understanding of the water balance is essential to explore water-saving techniques. In the 
context CRP project, experimental setups were conceived to monitor seasonal water consumption on the wheat crop 
irrigated by flood irrigation in Sidi Rahal station ( middle of morocco. The partitioning of evapotranspiration compounds 
shows that transpiration dominates the evaporation about 68 % for three days (22, 23 and 24 February 2012). In addition 
the wheat absorbs the soil water from 10 cm to 20 cm (90%) at this growing stage according to the multiple -source mass 
balance assessment. 
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I. INTRODUCTION 

With dramatic consequences, the drought, and meteoric characteristic, is one of a principal and structural data for Moroccan 
agriculture. The arid and semi arid areas take up 27% of Moroccan surface and 87% of arable land, whose 60% is covered by 
cereal culture [1], The cereal production remains insufficient to meet the country interior needs in this food product. To 
increase productivity by the improvement of the farming techniques is not only imperative, but it is the only solution to be 
considered. Enormous efforts were made to increase the yield. But, the techniques used are limited by drought which affects 
this crop during all its vegetative cycle. The irrigation of cereals remains the only possible way to improve the production. 
According to the FAO estimation, until 40% of food products produced in world are cultivated under irrigation, however a 
great quantities of water used for this purpose are lost by the escapes in the irrigation systems. In addition, the irrational 
irrigation practices are moreover one of the principal soil salinity causes. Approximately, 1/10 of the irrigated surfaces in the 
world were degraded by salt. Also the climatic changes make more and more areas in world exposed to the drought and 
desertification risk. The improved irrigation practices will contribute to preserve water and to protect the vulnerable soil. One 
widely used approach is conventional deficit irrigation (DI) but it requires crop specific information for its effective use [2], 

In the arid and semi arid areas, which suffers from a water shortage due to the scarcity of rainfall and to the increasing 
demands of water under the demography pressure and agricultural activity effect, the evapotranspiration constitutes the most 
important factor of water loss, whose, its determination is capital for a well control for water resources management. 

Contrary to classic methods traditionally used, which remain insufficient (micro -lysimetry, sap flow) to determine correctly 
the evapotranspiration partition, the isotopic geochemical studies provides important information’s to conclude this partition 
and to understand the extraction water processes by roots. Indeed, the heavy stables isotopes content increases in water soil 
by soil evaporation. On the other hand, the water extraction by the roots is without effect on this concentration. Water from 
growing-season precipitation is rapidly lost from the rooting zone by the transpiration or soil evaporation depending, in part, 
on the size of the precipitation event [3] and structural and physiological characteristics of the vegetation [4] the ratio of 
transpiration (T) to evapotranspiration (ET) is a synthetic parameter that integrates ecophysiological and micro 
environmental controls on total ecosystem water exchange [5]. 
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Indeed, in an isotopic steady-state condition of leaf water, transpiration introduces into the atmosphere a vapour whose 
isotopic signature is identical to that of root water [6], In a 5 2 H-5 18 0 diagram the signature of water vapour originally from 
transpiration belongs to the local meteoric water isotopic composition. The evaporation causes a little modification of 
isotopic composition of rainfall seepage in the surface layers of ground [7, 8], Moreover, the roots system, which is often 
widely developed in these surface soil layers, allow, by evaporation the plant alimentation in heavy isotope rich water. At 
last, the root extraction changes strongly the water distribution in different soil horizons and then the availability towards the 
vegetable covers. There are obviously a lot of consequences on the agronomic and hydrologic levels. The soil relative 
humidity, the hydraulic potential, the hydraulic conductivity, the root structure, the chemistry of soil solutions, and the 
evaporation - transpiration ratio are a parameters whose the determination is essential to understand the water transfer in soil- 
plant-atmosphere continuum. 

In this regard, a better understanding of the water balance is essential for exploring water-saving techniques and to avoid the 
contamination of ground water. The most important components of water balance in semi -arid areas are the 
evapotranspiration and the deep percolation. Effective methods, such as lysimetric method, sap flow measurement 
techniques, and micrometeorological techniques are used to measure or estimate ET. But, there are several limitations in 
using these methods. Stable isotopic tracer methods offer a new opportunity to study the components of ET at the field-scale, 
from the leaf level to ecosystem, and can partition the ET from different compartments of the ecosystem incorporating 
measurement of water vapour. 


II. Materials and Methods 

2.1 Study site 

The study is conducted in experimental station, irrigated with flood irrigation system; the experiment took place in one of the 
irrigated areas, which has been managed since 1999 by a regional public agency (Office Regional de Mise en Valeur 
Agricole du Haouz (ORMVAH)). The fields are generally sown between November 15 and January 15, and the harvest 
occurs after 5-6 months, in May or June. The station is located approximately 30 km southwest of Marrakech city, Morocco. 
This area has a semi-arid Mediterranean climate, characterized by low and irregular rainfall with an annual average of about 
240 mm against a higher reference evapotranspiration (ET 0 =1600 mm/year). The soils have high sand and low clay contents 
(18% clay, 32% silt, and 50% sand). 

2.2 Meteorological data 

The study site was equipped with a set of standard meteorological instruments to measure wind speed and direction (model 
Wp200, R.M. Young Co., Traverse City, MI, USA) and air temperature and humidity (model HMP45AC, Vaisala Oyj, 
Helsinki, Finland) at four heights. Net radiation over vegetation and soil was measured using net radiometers (a model 
CNR1, Kipp and Zonen, Delft, The Netherlands and the Q7 net radiometer (REBS Inc., WA, USA)). Soil heat flux was 
measured using soil heat flux plates (Hukseflux). Water content reflectometers (CS616, Campbell Scientific Ltd.) were 
installed at depths of 5, 10, 20, 30, 40, 60 and 80 cm in order to measure the soil humidity profile. Measurements were taken 
at 1 Hz, and averages stored at 30-min intervals on CR23X data loggers (Campbell Scientific Ltd.). 

2.3 Stable Isotopes Measurements 

2.3.1 Soil water, plant water and vapor sampling 

Using a hand-auger, soil was sampled from the surface to 10 cm. Sampled branches of orange tree were 0. 5-1.0 cm in 
diameter, 1-2 cm in length and from each of them the bark was removed. Every plant sample was composed of 2-3 stems 
from different individuals. Soil and plant samples were placed into screw-cap glass vials (5 ml) and sealed with Parafilm, 
then stored at about 2°C. Soil and plant water was extracted by cryogenic vacuum distillation [10], 

Water vapor was collected from 4 heights at the same time. During the collection period mentioned above, sampling was 
started at 10:00, 11:00, 13:00, 14:00 and 15:00 h. For each group, vapour was collected during lhour with a flow rate of 250 
ml min" 1 using a vacuum pump. The air was circulated through a set of 45 cm long glass traps [11] which were immersed in a 
mixture of ethanol and liquid nitrogen (about -80°C). Traps were made of 9 mm diameter Pyrex glass attached to 6-9 mm 
diameter Cajon Ultra-Torr adapters which framed in 9 mm diameter Swagelok Union Tee. After sampling the traps were 
sealed with Parafilm and stored at about 2°C. 
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Near the vapor sampling inlets, probes of model HMP45AC, Vaisala Oyj, Helsinki, Finland for measuring the air 
temperature (T a , in Kelvin) and relative humidity (/;,) every 5 min. Using T a , h and atmospheric pressure (Pa, in /zPa), water 
vapor concentration (mmolmol 1 ) was calculated by [12]: Eq 1 

H 2 0 = 10/7[P fI exp(13.3 185f -1.9760 2 -0.6445f 3 -0.1299 t 4 )]/P a 0 ) 

Where P a is standard atmosphere pressure (about 1013.25 hP a ) and t = 1 (373.15/T a ). 

2.3.2 Stable isotope and data analysis 

In the laboratory Soil and plant water was extracted by cryogenic vacuum distillation [10]. The water samples were 
isotopically analyzed at National center of sciences and nuclear techniques by laser spectrometer DLT-100 (+1 standard 
deviation). The standard deviation for repeated analysis of laboratory standards was 0.2%o for l8 0 and l%o for D. 
Concentrations of these isotopes are expressed as deviation from an international standard (V -SMOW) and using S notation 
in per mil (%o): Eq 2 

8% o = [(Rs / Rst) - 1] * 1 000 (2) 

Where Rs and Pst are the molar ratio of the heavy to light isotopes in the sample and the standard, respectively. 

2.3.3 Theoretical overview 

The isotopic ratio of the atmospheric water vapor at a certain altitude can be described using Eq. (3) by considering mixing 
of evapotranspired water vapor and free atmospheric water vapor [13, 14]. This relationship is linear, and when used with 
water vapor the v-intercept reflects the source isotopic composition of the evapotranspiration flux: 

8 ,, = C (5 — S FT ) f S FT 

ebl a ' a ET ' ^ ET , o\ 

^ebl 


Where (5 eb , is the isotopic composition of vapor collected from the ecosystem boundary layer, C a the atmospheric vapor 
concentration, C e bi the vapor concentration in the ecosystem boundary layer, d d is the isotopic composition of the atmospheric 
background and d ET indicates the isotopic composition of the evapotranspiration flux. The Keeling plot approach is based on 
the assumption that the atmospheric concentration of vapor in an ecosystem combines the inputs of two major sources: the 
background vapor from the atmosphere and vapor added by the sources in the ecosystem. It is further assumed that the only 
loss of water vapor from the ecosystem is by turbulent mixing with the background atmosphere. 

The isotopic ratio of evaporated water vapor from the soil surface is described below by considering the fractionation process 
[15] Eq 4: 


= 


« * S surf ~ h8 aun - s ea - (1 - h)s 


atm eq 


(l-h) + (l -K)E k b \ 


1000 


(4) 


5 E : is the isotopic composition of soil evaporation flux 

a* : is the temperature dependent equilibrium fractionation factor 

s k : is the kinetic fractionation factor 

h : is the relative humidity normalized to the temperature at the evaporation surface in the soil 
5 atm : is the isotopic composition of atmospheric vapour 

^surface • is the isotopic composition of water at the evaporation surface in the soil. 

The best equations describing a* are those provided by Majoube [16] Eq 5 

18 Oa* = [1 . 1 3 7(1 0 6 T 2 ) - 0.4 1 56(1 0 3 T) - 2.0667] / 1 000^ + 1 
Da* = [24.844(1 0 6 r 2 ) - 76.248(1 0 3 T) - 52.612]/1000+ 1 


(5) 
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Where, T is soil temperature recorded at 5 cm depth in degrees Kelvin. 

Ek is estimated using the diffusivity ratios of 1.0251 for H 2 0: HDO and 1.0281 for H 2 0:H 2 18 0 [17]. 
The contribution of transpiration to evapotranspiration is estimated by Eq 6 [18]: 


s t 


S T 



(6) 


III. Results and discussion 

3.1 Evolution of climatic conditions 

The temporal evolution of ET 0 is typically of a semi arid continental climate type. It is characterised by a high climatic 
demand. The lowest values of ET 0 (0.99mm/day) occurred during the winter and highest values occurred in the summer (9,89 
mm/day) in this figure the evolution of ET 0 is similar to T mean of air. Precipitation temporal patterns over the growing 
season of citrus trees were characterized by low and irregular rainfall events, with a total precipitation amount of about 295 
mm. The amount and timing of irrigations applied by the farmer are presented also in this figure (1) 



3.2 Partitioning Evapotranspiration components 

3.2.1 The stable isotopic composition of vapor water, soil water and irrigation water 

The following figure (2) chows that the continent of all samples for the three days of sampling (water vapor, soil water, 
irrigation water) are situated around the LMWL. The regression line of the all samples intersects the LMWL at point that 
present the origin of all the samples. The vapour profiles (figure 3) represent an isotopic heterogeneity during the three days, 
it gradually impoverishment -71.96%o to -87.01%o. However the vertical variation of 5D shows a homogenization in 
vegetation since 14:00 pm and it come on totality homogenization around 16:00 pm. That shows the high contribution of 
vegetation in hydrology cycle. 
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Figure 2: S 18 O versus SD in atmospheric water vapour, irrigation and soil water after observed at 

WHEAT STUDY CASE (A 22 st, B 23 th and C 24 th February) 



Figure 3 : Evolution of isotopic composition (6D)of atmospheric vapor at different levels in and out 

VEGETATION (B 23 th February 2012) 
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3.2.2 keeling plot analysis 

A significant regression lines were found in keeling plot pictured by 8D, although those plotted by 5 ls O were not significant. 
The table 1 shows the different slopes and the intercepts of keeling plots. All the intercepts of Keeling plots (figure 4) were 
close to symbols which reflected isotopic values of plant transpiration relative to that reflected isotopic compositions of soil 
evaporation. This meant plant transpiration contributes more to ET than soil evaporation. Considering wheat transpiration as 
one source and the soil evaporation as another one, the fractional contributions of plant transpiration to total ET (T/ET) were 
73 %, 59 % and 74 % for <5D on 22 st, 23 th and 24 th February, respectively. This may be explained by height canopy cover, 
about 77 % at this growing stage (maturation). The maturation period had larger E/ET ratios because of the smaller LAI and 
more bare area with low transpiration [19]. 

Yucui Zhang [20] combined isotopic and micrometeorologic approach to investigate the responses of transpiration and soil 
evaporation to an irrigation event in an winter wheat field, and the results show that transpiration was 4.91 mm, or 83%of 
total ET on DOY138 and 1.02 mm, or 60%of ET, on DOY149.The higher percentage transpiration during the filling stage is 
expected because this stage corresponds to a higher LAI and increasing biomass. 

Liu et al. [19] found that transpiration took up 70.3 and 69.7% of the total evapotranspiration for irrigated winter wheat and 
summer corn field in the growing season, respectively, at Luancheng Station in the North China Plain, based on the 
measurement of large-scale weighing lysimeter and two micro-lysimeters. These research conclusions about T/ET in 
irrigated field are similar to our results in this study, so the method of Linear mixing with atmospheric vapor (Keeling plots) 
to partition the evaporation and transpiration in study area is credible. 

However, for the researches not related to irrigation, the partitioning results have some difference. The T/ET ratio is 
relatively low for plants under semi-arid or arid climate [21]; [22], and relatively high under humid climate, which has 
similar moisture conditions as after irrigation. Hsieh [23] studied the T/ET in soils along an arid to humid transect in Hawaii 
by oxygen-18, and the T/ET ratio increased from 14%to 71%. 





Figure 4 : Relationship between SD and the inverse of the air absolute humidity at different levels 

ABOVE THE GROUND (A 22 st, B 23 th and C 24 th February) 
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Table 1 

Slope and intercept of the regression lines between AD values of water vapor collected at different 

HEIGHTS AND THE INVERSE OF THE CORRESPONDING VAPOR CONCENTRATION. THE INTERCEPT INDICATES THE ISOTOPIC 
VALUES OF EVAPOTRANSPIRATION (SET). 



5s 

5a 

5T 

SET 

P 

R 2 

a 

ck 

5E 

FT% 

22/02/12 

-27,26 

-73,46 

-34,55 

-57,58 

0,0002* 

0,530 

1,08 

1,03 

-119,28 

0,73 

23/02/12 

-5,72 

-68,22 

-27,29 

-52,95 

0,0001* 

0,562 

1,08 

1,03 

-90,09 

0,59 

24/02/12 

-26,05 

-71,36 

-37,44 

-59,75 

0,002* 

0,39 

1,08 

1,03 

-122,06 

0,74 


*The significance level is 0.05. 


3.2.3 Root water uptake depths 

For the soil and stem water sampled at the same time, according to the mass balance of soil water and its isotopes, the 
proportions of each layer (fl-f4) can be determined by their isotopic signature and the isotopic signature of mixture (stem 
water). Eq 7 

5Xt = f 15X1+ £25X2 + f35X3+f45X4 (7) 

With: 1= f 1 + f2 + O +f4 

IsoSource [24] was used to evaluate the relative contribution of each soil layer to stem water. The fractional increment was 
set at 1%, and the uncertainty level was set at 0.2. 



SOURCE: 2040 


SOURCE: 40-100 





Figure 5 : Histograms obtained by the multi-source mass balance, which show the estimated ranges of 

PROPORTIONAL CONTRIBUTION OF WATER FROM EACH SOIL DEPTH TO TOTAL WATER UPTAKE OF WINTER WHEAT 

(GREEN LINE STEM WATER, BLUE LINE IRRIGATION WATER) 
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The multiple-source mass balance assessment shows that the wheat absorbs the soil water from 10 cm to 20 cm (90%) at this 
growing stage (figure 5A) 

Figure 5B show hydrogen isotope in the soil profile and corresponding crop stem. The isotopic profiles of soil water are 
determined by the evaporation effect and the antecedent precipitation and irrigation. Because soil was sampled when there 
was no rainfall several days before, the surface soil water was isotopically enriched due to evaporation. For other soil water, 
the isotopes were a mixture of evaporation effect and the isotopic signatures of antecedent precipitation and irrigation event. 

In direct inference approach, the intersection of stem water isotopes ’vertical line and the soil water isotopes ’profile is 
considered to be the main depth of root water uptake. According to this direct inference we can get: the main depths of root 
water uptake of winter wheat are between 10 and 20 cm in this growing stage. 

IV. Conclusion 

This study, shows that the Keeling plot technique give a good result for the partitioning of Evapotranspiration components 
with high R 2 . The direct comparison shows: the main depth water uptake of winter wheat is 10 cm in this growing stage. The 
multiple -source mass balance assessment give the same result that direct approach, this one shows that the wheat absorbs the 
soil water from 10 cm to 20 cm (90%) at this growing stage. 

The present project constitutes an appropriate complement of previous attempts to quantify evaporation in irrigated areas 
within selected sites including different type of irrigation practices and different kind of vegetation (culture) and. This 
conclusion is useful for agricultural water management and irrigation schedules. 
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